Metal-nonmetal transition and excitonic ground state in InAs/InSb quantum dots 
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Using atomistic pseudopotential and configuration-interaction many-body calculations, we predict 
a metal-nonmetal transition and an excitonic ground state in the InAs/InSb quantum dot (QD) 
system. For large dots, the conduction band minimum of the InAs dot lies below the valence band 
maximum of the InSb matrix. Due to quantum confinement, at a critical size calculated here for 
various shapes, the single-particle gap Eg becomes very small. Strong electron-hole correlation 
effects are induced by the spatial proximity of the electron and hole wavefunctions, and by the 
lack of strong (exciton unbinding) screening, afforded by the existence of fully discrete OD confined 
energy levels. These correlation effects overcome Eg, leading to the formation of a bi-excitonic 
ground state (two electrons in InAs and two holes in InSb) being energetically more favorable (by ~ 
15 meV) than the state without excitons. We discuss the excitonic phase transition on QD arrays 
in the low dot density limit. 

PACS numbers: 71.35.Lk, 73.21. La, 73.22.-f 



The formation of excitons in semiconductors and in- 
sulators usually requires energy, e.g. photons, for one 
has to excite carriers across the single-particle band-gap 
Eg. There is a special interest, however, in the possibil- 
ity of forming excitons exothermically, i.e. an "excitonic 
ground state" as envisioned by Mott l] and Keldysh et 
al. 0- Indeed, the electron- hole system exhibits a rich 
range of phases 0, 0| as a function of the carrier density 
and effective-mass ratio TOg/m/i, including various exci- 
tonic insulating states such as molecular solid, exciton 
liquid, Mott insulator, and also various metallic phases. 
The excitonic ground state is of fundamental interest 
in itself because excitons can be a better alternative to 
atoms for studying Bose-Einstein condensation on 
account of the lighter excitonic mass, thus higher con- 
densation temperature. It is natural to search for exci- 
tonic ground states in systems where Eg is small, yet the 
screening is weak enough so as to prevent unbinding of 
the exciton. The search in hulk solids has thus focused 
on indirect gap semiconductors and semimetals to reduce 
screening, but excitonic ground states have not been con- 
clusively observed so far in such systems. Ground state 
excitons were also searched in nana structures, specifically 
in spatially indirect quantum- wells 0, , where electrons 
and holes are confined in different spatial regions. In 
"type 11" systems such as GalnAs/InP or CdTe/CdSe, 
electrons are localized in the well, whereas holes are lo- 
calized on the barrier, so screening is weak, but Eg is 
finite. In contrast, in "type III" heterostructures such 
as 01 InAs/GaSb, the conduction band minium (CBM) 
of the InAs well is lower than the valence band maximum 
(VBM) of the GaSb barrier, so at certain a well thickness 
one can have Eg—>0 well as separation of electrons 

from holes. Thus, at this thickness, one could expect an 
excitonic ground state if the electron-hole correlation en- 
ergy will be large enough to stabilize the complex. Re- 
cent experiments show evidence for existing excitonic 
ground state in type III InAs/Alj;Gai_2;Sb quantum-well 



superlattices; however, the binding energy was small (es- 
timated at ~ 3-4 meV) and strong magnetic fields are 
needed to stabilize the system. This is because in a 2D 
periodic well or superlattice, extended states exist in the 
in-pane direction leading to rather weak excitonic bind- 
ing. 

The advantage of a type-Ill system can, however, be 
utilized without the disadvantage of 2D periodicity un- 
derlying a quantum-well, if one considers a type III OD 
quantum-dot (QD). The well-known InAs/Al2;Gai_a;Sb 
epitaxial system [llj is inappropriate here, since it is 
lattice-matched, so dots will not form in a strain-induced 
Stranski-Krastanov (SK) growth 12] . We propose here a 
new dot matrix system, InAs/InSb, which has a type-Ill 
band alignment 13], and a lattice mismatch, and is hence 
amenable to epitaxial SK growth. Using an empirical 
pseudopotential approach [lj|, we find that as the dot 
size is reduced, electron levels localized on InAs move up 
in energy, whereas hole states, localized at the interface, 
initially above the electron levels, move down in energy. 
The system reaches a degeneracy Eg ~0 of electron and 
hole single-particle levels at a critical size, predicted here 
for realistic dot shapes. As a result of the OD confinement 
of both electrons and holes, there is but a small number of 
fully discrete bound states, so screening is limited. Using 
many-body configuration- interaction (CI), we find that 
consequently at the critical size, the exciton is bound by 
as much as ~15 meV ^ Eg, thus forming an excitonic 
ground state. Our study characterizes theoretically the 
properties of the excitonic ground state in this system, 
offering experimentally testable predictions. 

Identifying the material system: The customary way of 
selecting a strained dot/matrix system [l^ is to require 
that the dot material has a smaller gap than the ma- 
trix material (e.g., InAs/GaAs or CdSe/ZnSe) so that all 
carriers can be confined inside the dot. In conventional 
semiconductors Tsj the requirement Eg°^ < £;™atrix 
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plies a > 0,™^-'''!'=^ where a is the bulk lattice constant, 
leading to compressive strain. Here, however, we seek 
Eg^O with spatially separated, yet confined carriers, so 
we might need a system where the CBM of the dot ma- 
terial is below the VBM of the matrix material. This 
could be met when the dot material has a larger band 
gap (thus, a smaller lattice constant) than the matrix 
material, leading to tensile strain. SK growth under ten- 
sile strain was demonstrated before for PbTe/PbSefl^. 
The InAs/InSb system 0| has a lattice mismatch of ^ 
7%, and the CBM of InAs is below the VBM of InSb by 
0.17 eV for the unstrained system. This guarantees that 
the electrons will be confined on InAs, and holes on InSb. 
We will see below that one can further chose a strained 
dot shape that will bring the holes to the interface, thus 
providing proximity, yet spatial separation of electrons 
and holes. 

Strain and confinement: We have tested various dot 
shapes and sizes. Here we show three lens-shaped InAs 
dots of base/height =104/26, 143/36 and 156/39 A la- 
beled as Dl, D2 and D3, respectively, as well as spher- 
ical dots, all embedded in an InSb matrix. The su- 
percell that includes the dot and the matrix contains 
up to about 6x10^ atoms. The strain is relaxed by 
minimizing the strain-energy with respect to the dis- 
placement of all atoms, representing the energy as a 
sum of bond-bending and bond-stretching force con- 
stants (the valence force field method (VFF) 0). De- 
noting by eoif3, the a/S-component of the strain tensor. 
Fig. ^a) shows the strain profile along the (001) direc- 
tion across the dot D2, demonstrating that the isotropic 
strain I=exx+<^yy+^zz is positive (tensile), almost a con- 
stant inside the dots and decays to zero very rapidly 
in the barrier, whereas the biaxial component B — 

1/2 

[{exx - £yvf + {e-yy - ^zzf + {^zz " ^xxf\ IS nouzcro 
and nonconstant inside the dot, but decays slowly to 
zero outside. To observe how this strain modifies the 
electronic properties, we show in Fig. ^b) the strain- 
modified confining potentials obtained by inserting the 
strain of Fig. QJa) into the Pikus-Bir equations [l^. The 
isotropic strain lowers both the CBM and the VBM of 
InAs, but the lowering of CBM is greater, leading to a 
large reduction of the gap from the unstrained bulk value 
of 0.42 eV to 0.01 eV. Very importantly, because of the 
curved shape of the strained dot leading to a large biax- 
ial strain at the interface, the heavy-hole potential has 
a maximum at the InAs/InSb interface. We thus expect 
(and find below) that hole states will localize at the inter- 
face, rather than be extend throughout the barrier. This 
will enhance electron-hole binding. 

Method of calculation: We relax the atomic posi- 
tion, {R„_q} via the VFF method and construct the to- 
tal pseudopotential of the system V{r) by superposing 
the local, screened atomic pseudopotential Va{r) of all 
(dot-|-matrix) atoms and the nonlocal spin-orbit poten- 





(a) Strain profiles along (001) direction 
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(b) Strained band-offset along (001) direction 
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FIG. 1: (a) The isotropic and biaxial strain profiles along the 
(001) direction; (b) Strain-modified band-offset along (001) 
direction calculated from Pikus-Bir model, el, hh, Ih and so 
denote the band characters of the confining potential. 



tials Vso: V{r) = Vso + Y.n,a ^"(i" " R-n^a)- We use the 
Linear Combination of Bulk Bands (LCBB) method 19|, 
where the Hamiltonian — 1/2V^ + V{y) is diagonalized 
in a basis e,A(k)} of Bloch orbitals of band index 
m and wave vector k of material A (= InAs, InSb, 
GaAs), strained uniformly to strain e. This LCBB ap- 
proach Il 9( | p roduces accurate results for many nanos- 
tructures [20j , and greatly surpasses in accuracy the k • p 
method which limits the basis to to=VBM-|-CBM at 
k=0 only. The pseudopotentials of Ref. are used for 
InAs/InSb with minor modifications Many-body ef- 
fects are included via the configuration interaction (CI) 
method'i^l by expanding the total wavefunction in Slater 
determinants for single and bi-excitons formed from all 
of the confined single-particle electron and hole states. 
The Coulomb and exchange integrals arc computed nu- 
merically from the pseudopotential single-particle states, 
using the microscopic position-dependent dielectric con- 
stant '25| . Our CI approach is very similar to the Bethe- 
Salpeter equations |23|, except that in the latter case all 
the exchange integrals are unscreened, whereas since the 
exchange potential has a long-range component |2^ . we 
do allow its screening. 

Single-particle states: Diagonalization of the atom- 
istic single-particle pseudopotential Hamiltonian gives 
the confined single-particle electron and hole energy lev- 
els shown in Fig. [3 The corresponding wavefunctions 
for D2 are shown in Fig. |3| We see that the four InAs- 
confined electron states el, e2, e3, e4 move down in en- 
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ergy as the dot size increases; whereas the interfaciahy 
confined (see. Fig. ^ hole states hi and h2 move up in 
energy as the size increases. For the smallest dot Dl, 
there is a finite, single-particle gap of 96 meV between 
the first confined electron and hole states. For dot D2, 
the gap reduces to ~6 meV. For the large dot D3, the 
lowest InAs-confined electron level lies below the VBM 
of the InSb barrier and hybridizes with it. The single- 
particle gap is thus negative, signaling a charge transfer 
from InSb to In As. 

Many-body states: The energies were carefully con- 
verged by increasing the basis in the many-body expan- 
sion. For an uncorrelated electron-hole pair, the energy 
is Eg, whereas correlation could reduce the energy of 
the monoexciton Ex by Ax = Eg — Ex- Bi-excitons 
Exx could be bound with respect to two excitons by 
Axx = 2i?x ~ Exx- The CI energies are shown in 
Fig. H The results show that for dot Dl, Ex=79 meV 
and £'xx = 161 meV both positive i.e, correlation energies 
do not overcome the single-particle gap. However, for dot 
D2 we find that correlation reduces the energy of the sin- 
gle electron-hole pair from 6 meV to Ex — —9 meV, 
whereas for two electron- hole pairs, the reduction is from 
12 meV to Exx = —15 meV. Thus, both mono- and bi- 
excitons are now stable, the bi-exciton being the ground 
state, lower in energy than the "no exciton" case of a 
fully occupied valence bands and fully empty conduction 
bands state. To see whether the shape of the dot affects 
these conclusions, we performed similar calculations for 
spherical InAs/InSb QDs. The critical diameter for the 
metal-nonmetal transition is now around 2i?=64 A. The 
calculated exciton and bi-exciton energies are shown in 
Fig. EJb) for two diameters (52 A and 64 A). For the 
2R—64: A dot, we find an exciton energy of Ex = —22 
meV, and a bi-exciton energy of Exx — —39 meV. Again 
the bi-exciton is the ground state. 

Properties of excitonic ground state: For lens-shaped 




FIG. 3: COLOR: Wavefunction of the InAs-confined electron 
states (el-e4) and the first two hole states (hl-h2). The trans- 
parent lenses indicate the positions of In As dots. The isosur- 
face enclose 50 % of the state density except for e4 which 
is only weakly confined. For e4, the isosurface enclose only 
about 10 % of the state density. The contour plot are slices 
of the density taken from choosen planes. 



dots, the monoexciton binding energy Ax are 17 and 15 
meV for dot Dl and D2 respectively, whereas the neg- 
ative bi-exciton binding energy of Axx = — 3 meV for 
both Dl and D2, indicating an undbond bi-exciton. For 
the spherical dot with 2i?=64 A , Ax = 19 meV and 
Axx — —5 meV. The bi-exciton is also unbound. The 
bi-excitonic ground state is similar to a helium atom in 
the sense that there are two positive and two negative 
Coulomb bound charges. These states can be identified 
as "excitonic molecule" states ^ . The calculated perma- 
nent electric dipole moment p = / x.p{:x.)d^x for the exci- 
ton and bi-exciton of the lens-shaped dot D2 are aligned 
along (001) direction (growth direction), and amount to 
Pz=3.58 eA for the exciton and ^^=7.16 eA for the bi- 
exciton. These dipoles could be measured experimentally 
and since the exciton binding energy is so high (~ 15 
meV), it can be studied at rather high temperatures. 

Given these excitonic dipoles, it is interesting to study 
the collective behavior of excitons in the QD arrays due 
to dipolar interactions. Recent developments in epitaxial 
growth made it now possible to grow arrays of QDs . 
At a low dot density, where the QDs are far apart, the 
excitons behave like in the single dot while couple only 
weakly via dipole-dipole interactions. For simplicity, we 
ignore the fine structure |20| of the excitons and assume 
there can be only one exciton per dot at most (the results 
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(b) Exciton and bi-exciton total energies 
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FIG. 4: Exciton and Bi-exciton energies of (a) lens-shaped, and (b) spherical, InAs dots embedded in InSb matrix. 



can be easily generalized to the case allowing two excitons 
per dot). The excitonic phase transition can be described 
using an Ising-like Hamiltonian, 

where, is the number of excitons on dot i, and its 
value can be either or 1. Ex(j]) is the formation en- 
ergy of an exciton on a single dot, which is a function 
of T] being dot size, density etc. The dipole-dipole in- 
teraction Jij between the exciton on site i and site j 
is a function of both dot density and array structures. 
We consider here the case where the average interaction 
(J) — i/2NJ2ijJij{v) < (ferromagnetic-like), where 
N is the total number of QDs in the array. At T=0 K, 
when E{r]) + (J) <0, all QDs are occupied by excitons, 
thus the ground state is an excitonic state, accompanied 
by a macroscopic dipole moment. When E{r]) + (J) >0, 
no exciton will form in the QDs lattice, resulting in the 
normal insulting state. Particular interesting is the quan- 
tum critical point, where E{r]) + (J)=0. Here the fully 
occupied excitonic state and the fully empty state have 
the same energy, and the state could be either excitonic 
or normal, which can be switched by applying an elec- 
tric field. However, at a critical temperature Tc, there 
are large fluctuations between these two states leading 
to large fluctuations of dipole moment. Furthermore, the 
size distribution and disorder of the QDs lattice may lead 
to more complicated "Bose glass" phase . 

At medium and high densities, where the QDs have 
stronger coupling, there could be numerous phases as a 
function of carrier density and mass ratios, as proposed 
in Ref. 0. In these regions, the picture of (bi-)exciton on 
a single QD is not valid anymore. However, the electron- 
hole density and mass ratio can be tuned via the pa- 
rameters of dots arrays, such as the size and confining 
potential of a single dot, the density of dots, and even 



the structure of the dots arrays, thus make them very 
promising for studying the electron-hole phase diagram. 

To conclude, our pseudopotential many-body calcula- 
tions predict a valence-to-conduction single-particle en- 
ergy level cross-over and a spontaneous formation of bi- 
excitonic ground state at the critical size of lens-shaped 
and spherical InAs/InSb dot system. There are two rea- 
sons that make this type-Ill QD system a very promising 
candidate for studying excitonic ground state. First, the 
electron-hole energy gap is tunable by the QDs size alone, 
and second, the excitons have much larger binding ener- 
gies than in the bulk, or QW due to the OD confinement 
for both electrons and holes. Experimental studies on 
such a proposed excitonic ground-state are called for. 
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